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Abstract 

We calculate the string tension and part of the mass spectrum of SU(4) and SU(6) gauge 
theories in 2+1 dimensions using lattice techniques. We combine these new results with 
older results for N c = 2, ...,5 so as to obtain more accurate extrapolations to iV c = oo. The 
qualitative conclusions of the earlier work are unchanged: SU(iV c ) theories in 2+1 dimensions 
are linearly confining as N c — > oo; the limit is achieved by keeping g 2 N c fixed; SU(3), and even 
SU(2), are 'close' to SU(oo). We obtain more convincing evidence than before that the leading 
large- A^ c correction is 0(1/ N^). We look for the multiplication of states that one expects in 
simple flux loop models of glueballs, but find no evidence for this. 



1 Introduction 



As a by-product of our recent work on fc-strings in D=2+l SU(4) and SU(6) gauge theories 
U we also calculated the lightest masses with J PC = ++ ,0 ,2 ++ ,2 quantum numbers. 
In this note we combine these new results with the previous work of one of the authors so 
as to obtain more accurate information about the N c — oo limit and the N c dependence of 
SU(iV c ) gauge theories in D=2+l. 

This note is intended as an addendum and we refer the reader to @ as well as our related 
work in D=3+l |1|, |3| for a discussion of the (well-known) physics motivation M as well as 
the (mostly standard) lattice techniques. 

In the next section we discuss the way our analysis differs from that of |§ and we perform 
a partial re-analysis of that earlier work so as to allow us to include it in a combined analysis 
with the new calculations. We then turn to an analysis of the string tension, a, which is the 
physical quantity we calculate most accurately. We confirm that the N c = oo limit is achieved 
by keeping constant the 't Hooft coupling, g 2 N c , and we are able to determine quite accurately 
that the leading corrections are 0(1 /N 2 ). We then analyse the glueball masses and confirm 
that they also show little variation with N c . Finally we describe some calculations designed 
to reveal the increasing multiplication of states with increasing N c , which is expected M in 
models where glueballs are composed of closed loops of flux. We find no evidence for this in 
our exploratory study. 



2 Some technical details 

Our lattices are cubic with periodic boundary conditions. We use the standard plaquette 
action 

S = /?£{! - ^ReTrU p } (1) 

where the ordered product of the SU(iV c ) matrices around the boundary of the plaquette p is 
represented by U p . In the continuum limit this becomes the usual Yang-Mills action with 

where we recall that g 2 has dimensions of mass in D=2+l, so that the dimensionless bare 
coupling, being a coupling on the scale a, is just ag 2 . The simulations are performed with 
a combination of heat bath and over-relaxation updates, as described in 0. Our SU(4) and 
SU(6) updates involve 6 and 15 SU(2) subgroups respectively. Masses are calculated from the 
exponential fall-off of appropriate correlators at large enough t, again as described in 0. In 
Table 1 we list the values of (3 at which we perform calculations, together with the lattice size, 
L, and the value of the average plaquette. In Tables 2 and 3 we list the masses we obtain for 
SU(4) and SU(6) respectively. 

Our present calculations differ from || in some technical details. Firstly, the mass cal- 
culations have been performed for fewer quantum numbers and for a more limited range of 
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operators. (The primary purpose of these calculations was to obtain the /c-string tensions dis- 
cussed in |IJ.) Since the mass calculations involve a variational selection of the 'best' glueball 
wavefunction from the basis of operators used, there is a danger that this might introduce a 
problem in combining the present results with those of @. However we shall see that the new 
and old SU(4) calculations agree with each other within errors which tells us that there is in 
fact no reason to worry. 

The second difference is that we now do exponential fits to the correlation functions (as 
described in ||) rather than relying on effective masses as in [El. Again, the fact that the 
new and old SU(4) calculations agree with each other within errors reassures us that this 
introduces no significant relative bias. 

Finally we use standard Gaussian error estimates rather than the somewhat idiosyncratic 
algorithm described and used in [EJ . In order to be able to perform a unified statistical analysis 
of the old and new calculations, we have redone all the (relevant) continuum extrapolations 
in 0, obtaining the values shown in Table 4. We note that the new errors are typically about 
30% smaller than those presented in |EJ. 



3 The string tension 

Our first step is to calculate the mass, mp(L), of a flux tube that winds around our L x L 
spatial torus, as described in |2j, |]. We then extract the string tension, a, using 

amp(L) = a 2 cxL . (3) 

6L 

This relation assumes that we have linear confinement and that the leading correction is that 
of a simple periodic bosonic string ||. In the case of SU(2) there is now rather convincing 
evidence fl| for both these assumptions. As an example of what one can do for larger groups, 
such as SU(4), we plot amp(L) as a function of L at (3 = 28 in Fig.l. (The masses on the 
smaller volumes are listed in Table 11 of [Q.) We see the approximate linear rise confirming 
the persistence of linear confinement for N c = 4. (Note that the largest string length is 
16a ~ 4/y/a which should be long enough to realistically address this issue.) We also see that 
the bosonic string correction fills in most of the deviation from linearity at larger L. We shall 
therefore use eqn(^) to extract the string tension from the masses of flux loops that are longer 
than ~ 3/y/a, where one typically finds |TJ that eqn([3]) works well. 

Having calculated the string tension as a function of (3, we follow [Q] and extrapolate to 
the continuum limit using 

fray/a = c + ^- (4) 
Pi 

where (3j is the mean field improved coupling 

/ 3 I = f3x(±-TrU p ). (5) 

We obtain the continuum string tension from the fitted value of cq = 2N c y/a/g 2 . The values 
of yjo 1 ' g 2 thus obtained are listed in Table ||. 
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In Fig.|2] we plot the calculated values of \Ja / g 2 against N c . We immediately observe that 
the variation approaches a linear form for larger N c 

^ocN c (6) 

r 

and indeed is nearly linear even down to N c = 2. This confirms the usual diagrammatic 
expectation [|J that a smooth large-iV c limit is obtained by keeping constant the 't Hooft 
coupling A = g 2 N c . 

We can also ask what is the power of the leading correction. To do so we fit our string 
tensions with the functional form 

c o + TfZ- (7) 
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In FigJJ we show how the goodness of fit varies with the power a. From this we can infer that 

a = 1.97 ±0.21 (8) 

if we choose errors corresponding to a 20% confidence level. If we assume that the power 
should be an integer, then only one value is allowed: a = 2. This is in agreement with the 
usual diagrammatic expectations [|J. 

Fitting our calculated values with a = 2, we obtain 

> - 0,9755,34) - «£) (9) 
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This fit has a good confidence level, ~ 70%. Thus, despite our smaller errors and extended 
range in N c , the observation in |2| still stands: we can describe the string tension of SU(N C ) 
gauge theories, all the way down to SU(2), by that of the SU(oo) theory supplemented by 
the leading correction with a modest coefficient. 



4 The mass spectrum 

We list in Tables ^] and ^] our calculated masses. We know |7J that with the plaquette action 
the leading lattice correction to dimensionless ratios of physical quantities, such as mo/ ' \Jo is 
0(a 2 ). Thus we extrapolate to the continuum limit using 

m G {a) m G (a = 0) i 2 



y/a(a) yja{a = 0) 



+ caV (10) 



The results are shown in Table ^. All these extrapolations have a good x 2 apart from those 
indicated by an asterisk. We also show similar extrapolations for the masses calculated in 
but with the errors recalculated as described in Section 2. 

We first note that the old and new SU(4) calculations agree, within errors, for the lightest, 
most accurately determined states. This reassures us that the different procedures used in 
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these calculations do not lead to a relative systematic difference that would undermine a 
combined statistical analysis. 

The variation of our mass ratios with N c is weak which encourages us to try to fit with 
just the leading large- N c correction: 

m G , d x 

d o + T7t (11) 



9 2 N C iV c 2 

where do = \im.^ me / g 2 N c . The results of such fits are listed in Table ^ where we also 
produce corresponding values of mc/^cr, which have been obtained by using eqn(||). We note 
that we obtain acceptable fits in all cases, except in the case of the rather heavy 2 ++ where 
the low value in SU(6) makes it impossible to find a good fit. 

To illustrate all this we plot mg/ ' g 2 N c against 1/N 2 for the ++ and glueballs in Fig.|j. 
We observe that the dependence of our masses on N c is really very weak over the whole range 
ofAL 



5 Are glueballs loops of flux? 

If SU(3) glueballs are composed of closed loops of fundamental (k=l) flux then in SU(iV c > 4) 
there will be extra states composed of closed loops of k ^ 1 flux |J. In simple, but natural, 
string models of glueballs |5], || the masses of such states will be scaled up by a factor of 
approximately yak/a, where <j\. is the fc-string tension (so that /Jk=i = c). These string 
tensions have been recently calculated and satisfy Casimir scaling, <7fc/cr ~ k(N c — k)/ (N c — 1), 
to a good approximation Is there any sign of such extra states in our calculations? 

We recall that our glueball operators are composed of products of highly 'smeared' SU(iV c ) 
matrices around the boundaries of squares and rectangles @. When we take products of 
these matrices around the torus we obtain a basis of operators that typically provides a good 
wavefunctional for the winding fundamental flux tube and so we may regard our usual glueball 
operators as being closed loops of fundamental flux. To construct corresponding loops of k = 2 
strings we take the fundamental loop, I, twice around the boundary in the form Tr(Z.Z) or Tr/. 
We can project onto the antisymmetric k = 2 representation, k = 2A, as well. We look for 
extra states in the k = 2 and k = 2A bases of operators as compared to the k = 1 basis. We 
note that the lightest few states obtained in the usual k = 1 basis seem to continue well from 
N c = 2, 3 to larger N c and so seem to possess no extra states at larger N c . 

We calculate the spectrum of states separately using as a basis either the usual k — 1 
operators, or the k = 2 operators or the k = 2A operators. We do this for SU(6) where mixings 
and decays, which might serve to obscure the effects we are searching for, should be weak. 
In SU(6) any such new states should have masses that are larger by a factor ~ \/1.6 ~ 1.26. 
In Table [7] we show the ++ states obtained using these different bases of operators. Masses 
are obtained from exponential fits. For heavier states, indicated by a single asterisk, we take 
the effective mass from t — a to t — 2a. for the very heaviest states, indicated by a double 
asterisk, we take the effective mass from t — to t — a. What we see is that there is no clear 
sign of extra states in the k = 2A sector as compared to the k = 1 sector, in the expected mass 
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range. It is interesting to note that the one clear extra state in the k = 2 sector, with a mass of 
am ~ 1.20, fits in well with being just the scattering state composed of two ground state ++ 
glueballs in an S'-wave. (Large- N c counting will suppress its appearance in the k = 1 sector.) 
We have performed similar comparisons in the and 2 ++ sectors with similar conclusions. 

Although our results look quite negative, they should be regarded as preliminary. Our 
operators are based on closed loops whose sizes are not much larger than the transverse extent 
of the smeared matrix, and which might therefore not be a good basis for closed 'string-like' 
states. So calculations with larger loops need to be performed. Moreover the lighter states 
may be smaller and may be the least string-like. Thus one needs an accurate comparison for 
higher mass states, which is something we would not claim to have achieved in the present 
calculation. 



6 Conclusions 

In this paper we have confirmed the main conclusions of over a larger range of N c . We 
have also looked for an increase in the density of states of the kind that naturally arises 0, |S| 
in models where glueballs are composed of closed loops of flux. In what is admittedly a crude 
study, we found no sign of such extra states. 

One motivation for these large- N c calculations is to provide theorists who are attempting 
to develop (mainly) analytic methods, with something to compare against. Thus we find it 
interesting to note that the string tension derived as a first approximation in M 



" (12) 
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is in fact extremely close to all our calculated values. As for the glueball spectrum, we note the 
predictions in |IJ for the lightest glueball masses in the N c — * oo limit: m ++ = 4.10(13)y / a, 
m — /m ++ = 1.35(5) and m 2 ++/m ++ = 1.60(17). The first and last of these predictions are 
consistent with our values, while the second is only two standard deviations out. We also note 



that recent AdS/CFT calculations of glueball masses (see JTTJ] for a review) are compatible 



with our values. Finally one should add that insights into the significance of the detailed mass 



spectrum may be obtained through a comparison with related D=2+l field theories [12 
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SU(4) 


SU(6) 


P 


L 


plaq 


P 


L 


plaq 


60.0 


32 


0.914426(2) 


108.0 


24 


0.887953(2) 


45.0 


24 


0.884724(3) 


75.0 


16 


0.835392(6) 


33.0 


16 


0.840190(7) 


60.0 


12 


0.790321(8) 


28.0 


16 


0.809332(6) 


49.0 


10 


0.736809(12) 


28.0 


12 


0.809350(15) 


42.0 


8 


0.684495(20) 


21.0 


10 


0.737640(15) 








18.0 


8 


0.68593(4) 









Table 1: Average SU(4) and SU(6) plaquette values for the new calculations. 



state 


(3=18 


(3 = 21 


(3 = 28 


(3 = 28 


(3 = 33 


(3 = 45 


(3 = 60 




L=8 


L=10 


L=12 


L=16 


L=16 


L=24 


L=32 


/ 


1.50(2) 


1.223(11) 


0.715(4) 


0.9937(58) 


0.6622(29) 


0.4974(22) 


0.3571(14) 


0++ 


1.85(6) 


1.54(4) 


1.057(15) 


1.082(11) 


0.872(10) 


0.620(5) 


0.4592(27) 


++* 


2.5(3) 


2.29(15) 


1.53(3) 


1.56(4) 


1.289(23) 


0.940(8) 


0.682(4) 


0++** 




2.58(34) 


1.77(5) 


1.88(6) 


1.506(24) 


1.117(33) 


0.842(5) 


0~ 


2.7(5) 


2.11(14) 


1.61(3) 


1.53(5) 


1.306(22) 


0.895(20) 


0.6736(34) 


* 




2.16(26) 


1.97(9) 


1.96(14) 


1.65(4) 


1.233(16) 


0.852(14) 


2++ 




2.4(3) 


1.5(2) 


1.81(5) 


1.502(26) 


1.051(12) 


0.773(6) 


2++* 






2.35(16) 


2.23(11) 


1.71(6) 


1.22(6) 


0.887(18) 


2~ 






2.03(11) 


2.21(22) 


1.81(5) 


1.228(13) 


0.921(10) 



Table 2: Masses of glueballs, labelled by J , and of the periodic flux loop, I, of length L : 
in SU(4) at the indicated values of (3 on lattice volumes L 3 . 
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state 


p = 42 


p = 49 


p = 60 


p = 75 


o 1 no 

p = 108 






t m 
L=1U 


T 1 O 

L=lz 


L=lo 


T O A 


/ 


1.453(14) 


1.194(9) 


0.8575(47) 


0.6825(31) 


0.4552(13) 


0++ 


1.754(54) 


1.491(31) 


1.161(16) 


0.8885(42) 


0.5925(26) 


o++* 


2.7(7) 


2.13(13) 


1.752(60) 


1.347(18) 


0.8769(57) 


Q++** 




2.58(34) 


1.97(9) 


1.605(27) 


1.1226(83) 


0~ 


2.50(27) 


2.20(13) 


1.707(43) 


1.284(17) 


0.861(8) 


* 




2.8(4) 


2.14(11) 


1.66(4) 


1.119(12) 


2++ 


2.6(4) 


2.6(4) 


1.92(7) 


1.478(23) 


0.953(16) 


2 ++* 






2.40(16) 


1.82(6) 


1.242(11) 


2— 






2.3(2) 


1.79(5) 


1.197(10) 



Table 3: Masses of glueballs, labelled by J PC , and of the periodic flux loop, I, of length L : 
in SU(6) at the indicated values of (3 on lattice volumes L 3 . 



m G /y/a 


state 


SU(2) 


SU(3) 


SU(4) 


SU(5) 


SU(4) 


SU(6) 


0++ 


4.716(21) 


4.330(24) 


4.239(34) 


4.180(39) 


4.235(25) 


4.196(27) 


++* 


6.78(7) 


6.485(55) 


6.383(77) 


6.22(8) 


6.376(45) 


6.20(7) 


Q++** 


8.07(10) 


8.21(10) 


8.12(13) 


7.87(18) 


7.93(7) 


8.22(12) 


0~ 




6.464(48) 


6.27(6) 


6.06(11) 


6.230(44) 


6.097(80) 


* 




8.14(8) 


7.84(13) 


7.85(15) 


8.20(15)* 


7.98(15) 


2++ 


7.81(6) 


7.12(7) 


7.14(8) 


7.15(12) 


7.17(8) 


6.67(18) 


2 ++* 






8.50(17) 


8.56(15) 


8.06(22) 


8.89(20) 


2— 




8.73(10) 


8.25(21) 


8.25(18) 


8.49(13) 


8.52(20) 



Table 4: Glueball masses in units of the string tension, in the continuum limit. Reanalysis 
of [Q] on left; new calculations on right. 
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group 


c 


Cl 


V*/9 2 


SU(4) 
SU(6) 


6.072(16) 
14.000(38) 


-8.19(45) 
-46.5(2.6) 


0.7590(20) 
1.1667(32) 


SU(2) 
SU(3) 
SU(4) 
SU(5) 


1.3405(31) 
3.3182(61) 
6.065(21) 
9.657(38) 


-0.417(24) 
-2.43(11) 
-7.74(73) 
-21.4(1.9) 


0.3351(8) 
0.5530(10) 
0.7581(26) 
0.9657(38) 



Table 5: Continuum extrapolations of fiia^/a — ► 2N c ^Ja/g 2 as in eqn(^). New calculations 
(top) and reanalysed old calculations (bottom). 



state 


Yim Nc _ >00 m/g 2 N c 


slope 


iV c > 




limjv^oo mj sja 


0++ 


0.8116(36) 


-0.090(28) 


2 


0.9 


4.108(20) 


++* 


1.227(9) 


-0.343(82) 


2 


0.8 


6.211(46) 


Q++** 


1.65(4) 


-2.2(7) 


4 


1.6 


8.35(20) 


o— 


1.176(14) 


0.14(20) 


3 


0.2 


5.953(71) 


* 


1.535(28) 


-0.35(35) 


3 


1.0 


7.77(14) 


2++ 


1.359(12) 


-0.22(8) 


2 


2.9 


6.88(6) 


2 ++* 


1.822(62) 


-3.9(1.3) 


4 


0.4 


9.22(32) 


2~ 


1.615(33) 


-0.10(42) 


3 


1.0 


8.18(17) 



Table 6: The large N c limit of the mass spectrum in units of g 2 N c ; with the slope of the 
linear fit when plotted against 1/N 2 . Also the range of colours fitted and the x 2 P er degree 
of freedom of the best fit. Last column uses the value of lim^^oo y/a/g 2 N c in eqn@. 
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am eff ; J FC = 0++ 


k — 1 ; 10 ops 


k = 2A ; 10 ops 


k = 2 ; 20 ops 


0.592(4) 


0.594(4) 


0.593(4) 


0.88(1) 


0.90(1) 


0.89(1) 


1.14(1) 


1.17(1) 


1.13(1) 






1.20(1) 


1.30(2) 


1.29(6) 


1.31(2) 


1.47(2)* 


1.42(2)* 


1.46(2)* 


1.53(2)* 




1.49(2)* 


1.78(4)* 


1.76(4)* 


1.72(4)* 






1.73(4)* 






1.65(3)* 






1.72(4)* 


2.9(2)** 


2.9(2)** 


2.3(1)* 



Table 7: Spectrum of scalar masses at (5 — 108 in SU(6) obtained using operators based on 
fundamental flux loops (k — 1) doubly charged loops (k = 2) and from the antisymmetric 
combination of the latter (k = 2 A). States are listed in the order of effective masses from 
t — to t — a and are matched to guide the eye. 
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ami(L) 



0.8 



0.6 



0.4 - 



0.2 - 




12 16 



L 



Figure 1: The mass of a periodic flux loop of length aL at (3 = 28 in SU(4). The line shows a 
fit of the form in eqn(|3[). 
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1.2 



1 - 



\J<y 

IF 



0.8 



0.6 - 



0.4 - 



0.2 




1 2 3 4 5 6 



Nr. 



Figure 2: The value of \jo j g 2 as a function of N c . The line shows the fit in eqn@. 
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Figure 3: The x 2 P er degree of freedom against the power, a, of the leading large- N c correction 
when fitting y/a/g 2 N c to eqn(|7]). 
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1/N t 



Figure 4: ++ («) and (x) glueball masses in units of g 2 N c , plotted against 1/N%. The 
best linear extrapolations to the N c = oo limit are shown. 
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